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Abstract. We show that in the large matrix Umit, the eigenvalues of the normal ma- 
trix model for matrices with spectrum inside a compact domain with a special class of 
potentials homogeneously fill the interior of a polynomial curve uniquely defined by the 
area of its interior domain and its exterior harmonic moments which are all given as 
parameters of the potential. 

Then we consider the orthogonal polynomials corresponding to this matrix model and 
show that, under certain assumptions, the density of the zeros of the highest relevant 
orthogonal polynomial in the large matrix limit is (up to some constant factor) given by 
the discontinuity of the Schwarz function of this polynomial curve. 
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Introduction 

In 1999, P. Wiegmann and A. Zabrodin realised [H, 0| that conformal maps from the 
exterior of the unit disc to the exterior of a simple closed analytic curve admit as functions 
of the exterior harmonic moments (t^)^]^ of the analytic curve the structure of an inte- 
grable hierarchy: the dispersionless Toda lattice hierarchy. (The Toda lattice hierarchy 
was introduced by K. Ueno and K. Takasaki in j^. A review of the dispersionless Toda 
lattice hierarchy can be found in [3].) 

As was shown earlier by L.-L. Chau and Y. Yu 0, Hj, the normal matrix model defined 
by the probability distribution 

(1) p^_^(M) = ^e-^*^^(*^), Zr,,v= [ P^,v(M)dM, 

on the space A/jv of all normal N x N matrices with the potential 




(2) V(M) = - I MM* - ^{tkM'' + tkM*^^ 

k=l 



gives rise to a solution of the Toda lattice hierarchy. In this connection between 
normal matrix models and conformal maps was picked up and the authors showed that 
the eigenvalues in the large matrix limit homogeneously fill the interior of the analytic 
curve with the exterior harmonic moments (tfc)fc^i and encircled area vrto- 

Nevertheless, all these results remain on the level of formal manipulations as in general 
already the integral in ([T]) diverges. 

The purpose of this work is now to proof these statements in a mathematically rigorous 
setting. To this end, we will first of all introduce a cut-off for the divergent integral by 
restricting ourselves to the space of those normal matrices whose spectrum lies in some 
compact domain. Furthermore, we will only consider polynomial potentials. Then, only 
finitely many exterior harmonic moments are different from zero and so the corresponding 
curve is polynomial. Since, at least for curves with small encircled area vrto, we have a 
bijective relation between a polynomial curve and its harmonic moments, we finally can 
prove that, provided to is small enough, in the large matrix limit the eigenvalues indeed 
homogeneously fill the interior of the polynomial curve uniquely determined by the exterior 
harmonic moments {tk)'^=i and the encircled area nto. 

Moreover, we find that the distribution of the zeros of the n-th orthogonal polynomial 
Pn,N associated to this matrix model is in the limit — x G [0, 1], ^ oo (up to some 
constant) given by the discontinuity of the Schwarz function of the polynomial curve with 
area irxtQ and exterior harmonic moments (tfc)fcLi) provided the zeros in this limit are 
confined to some one-dimensional tree-like graph. 

The paper is organised as follows: In chaper one we will introduce the considered matrix 



model and review a few standard methods (see e.g. [Ul]) to characterise the distribution 
of the eigenvalues. In the second chapter, we will show that polynomial curves with small 
area are uniquely determined by its exterior harmonic moments (for star-like domains this 
was already shown in [13] )• Then, after having introduced the concept of the Schwarz 
function (see e.g. we will end up with the connection between harmonic moments of a 

polynomial curve and its parametrisation, which we will see is given by the dispersionless 
Toda lattice hierarchy. Thereafter, we will consider the large matrix limit, in chapter 
three from a general point of view which is well know from potential theory 
and in chapter four, we will apply those results to potentials of the form with only 
finitely many no n- vanishing tk, leading us to the desired result that for small to the 
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eigenvalues homogeneously fill the interior of the polynomial curve determined by the 
exterior harmonic moments (tfc)^i and the area vrto- In the last chapter, we will analyse 
the behaviour of the orthogonal polynomials of these potentials in the continuum limit 
and will find that they obey the Toda lattice hierarchy. Finally, we will connect the 
limiting distribution of the zeros of the n-th orthogonal polynomial to the discontinuity 
of the Schwarz function of the corresponding polynomial curve and calculate it explicitly 
for the Gaussian and the cubic case. 

1. The model 

1.1. The eigenvalue distribution. We start with the probability distribution 

VNy{M) = J—e-'^^^viM)^ 2^^^ ^ f Vm,v{M) dM, 

on the space 

Un{D) = {M e C^^^ I [M,M*] = 0, a{M) C D} 

of all normal N x N matrices with spectrum in a compact subset D G G, where V(M) 
denotes some polynomial in M and M* such that V(M) = V(M)*. The measure dM on 
the variety A/Ar(-D) shall be the one induced by the flat metric on (D^><^. 

Proposition 1.1. On the regular part ofAfN{D), the measure dM factorises into 

N 

(3) dM=dU H \zi-z,\^l[dh„ 

l<i<j<N i=l 

where {zi}fL^ is the spectrum of M , M = U d\a.g{zi)f^^U* for U G U(iV)/U(l)^, and the 
measure 

(4) dU = H d\,,, d^i, = {U*),k d^Uk,, 

i<j 

on U(A^)/U(1)^ is induced by the measure on \J{N), which is given by the fiat metric on 

(^NxN 

Proof. Pulling back the tangent spaces of U(iV) to the Lie algebra u(iV), we uniformly get 
the metric g{X, Y) = tT{X*Y) on u{N). Since the tangent spaces of the orbits U ■ U(l)^, 
U G U(A^), are generated by the set d of diagonal matrices in u{N), the tangent spaces 
of U(A^)/U(1)^ are generated by u(A^)/0, and the metric g reduces to 

~g{X,Y) = 2J2^^,y^J, X,Yeu{N)/^. 

i<j 

With the upper triangular matrix elements as coordinates, we therefore get the measure 
(iD on U(A^)/U(1)^. 

Let now t ^ M{t) be a smooth curve in the regular part of A/7v(-D) with M(0) = M. 
Then we find a smooth curve 

t ^ {z{t), U{t)) eD^ X \J{N)/\J{lf with 
M{t) = U{t)D{t)U{ty, D{t) = diagU(t)), 
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which up to a permutation of the eigenvalues is uniquely determined. Again taking the 
pull back of the tangent space of U(A^), ii = U*U, the flat metric on C^^^ reads 



N N 

tr{M*M) = tT{D*D + 2u{D*u - uD*)D) = ^ z^Zi + ^ Uijii 

i=l *ii=l 



^ij I Zi Zj I 



Choosing ij and {uij)i^j as coordinates, we get the Riemannian volume form (jS]). □ 



N 



Since J\fi\i{D) is regular almost everywhere, we can integrate over the group U(A^) /U(l) 
leaving us with the probability distribution 

(5) PnA^) = ^e-^S"-^(-)|AU)r, Zr,,v = ( [ P^yU) d'^^zV , 

for the eigenvalues z = {zi)^i G , where A(z) = det{zi~'^)^j^i = Y[i<ji^j ~ ^i) denotes 
the van-der- Monde determinant and V{z) is the real- valued polynomial in z and z induced 
by V{z 1) = V{z) 1, 2 G Z). From now on, as long as it is clear which potential V we use, 
we will suppress it in the notation. 

Because the probability Pat vanishes if two eigenvalues are equal, we may consider it 
as probability distribution on the regular part 

= {zeD'' \z,^z^^i^j}. 

To describe the behaviour of the eigenvalues, we are going to introduce the notion of 
correlation functions. 

Definition 1.2. The k-point correlation function R^^\ 1 < k < N , is given by 



N\ r ^ 
(TV - A;)! Jdn-, .j-^-^^ 



So -^i?^^ = P/v, and j^R^^^ is the density of the eigenvalues. 

Definition 1.3. We denote by A]\r{n,X,a), < n < N , the probability of having exactly 
n eigenvalues in a disc Bx{a) of radius A > around some point a E D: 

(6) A^{n,\,a)=(J^ J PN{z)d?''z. 

B^(arx{D\B^(a))«-" 

Proposition 1.4. For Q<n<N,\>{], and a E D, we have 
Proof. Assume inductively that the equation 



(7) J f{z) d'-z= j^i-^r (i) J /U) d^"^ 
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which certainly is true for n = 0, holds for some A G D and any symmetric function / 
on for some n G Nq. Then, for any symmetric function / on 

f /U)d2("+i)^= f /U)d2("+i)^- f /U)d2("+i)^ 

Dx{D\A)" Ax(D\A)" 

k=0 ^ ^ ^"^fex£)n-fc + l J^k + ly,^n-k ^ 

n+1 / t\ r 

fc=0 ^ ^ + 1 

which proves the identity ([7]) for all n G Nq and any domain A G D. 
Therefore, 

A^in,X,a) = (^) (^r) [ ^^(^)d^'^^ 

V / fc=0 ^ ^ B;,(a)"+'=xD^-"-'= 

1.2. Orthogonal polynomials. One may formulate the problem of finding the correla- 
tion functions as the construction of a set of orthogonal polynomials. 

Definition 1.5. We call the polynomials Pn.N of degree n with leading coefficient one, 
n G No, defined with the inner product 

(8) {f.9)N= I 

JD 

by the relation 

{Pm,N,Pn,N)N = 5m,nhn,N for all m, n G Nq 

for some positive numbers hn^N, "H- G Nq, the orthogonal polynomials for the potential V 
with norms (/ira,Af)^o- -^^^ normalised functions 

qn,N = —^=Pn,N, n G Nq, 
V hn,N 

we call the orthonormal polynomials for the potential V. 

Proposition 1.6. Let {pn,N)'^=o be the orthogonal polynomials for the potential V with 
norms {hn,N)'^=o ^'^'^ 



Kn{w, z) = e-f (^M+^W) J2 ^"-^^^^^"-^(^^ 



n=0 ^"'^ 



Then 



(i) / KN{z,z)dz = N and / KN{w,v)K]y{v, z) d v = K]\j{w, z), 
Jd Jd 

(ii) Z^ = N\l[ hn,N and Rf{z) = det ^,)) J zgDK 



n=0 
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Proof. The first part is a direct consequence of tfie ortlionormality of tlie polynomials 
{—7^ — Pn,N)n=o- For ^lie second part, we write out the van-der- Monde determinants in 

equation ([5]), use the fact that the determinant is invariant under column operations, and 
get for z G 



PnU) 



Un=0 hn,N ^-N E£i V{z,) 

Zn 

nN-l , 
n=0 ^n,N ^-Ny,^_, V(z, 

Zn 

nN-l , 
n=0 '^n,N 



N 



^ sign(ra)JJ 

N /N-1 

t&Sn «=1 \n=0 

N 



Pa(i)-l,Ar(^i)PtT(i)-l,Af(^ra(i)) 



ha(i)~l,N 



Pn,N{Zi)Pn,N{Zr{i)] 



'■n,N 



Using the properties of part (i), we find 



Jd 



detiK^izi, Zj))ljl-^ I KN{zk,Zk)d^^Zk 



D 



k-1 



j = l creSk-i 



I i=i 



N\,Zij Zcr(i) J 



D 



KN{zj, Zk)K{zk, z^(^j)) d^Zk 



{N -k + l)det{KN{z,,Zj)fC^^ 



Integrating now recursively the expression for P/v, we get the relations of part (ii). □ 

Proposition 1.7. For n G N, the n-th orthogonal polynomial Pn,N for the potential V on 
the domain D is given by 



(9) 



Pn,N{zo) 



» n 

/ l[{zo-z,)P^Ny{z)d^''z, zoeC. 



Proof. By definition, the polynomial 



Pn,N[Zo) 



det (^^-^.r^e-^^(^.)d2^);_ 



which by construction is orthogonal to all monomials of degree less than n and has leading 
coefficient one, is the n-th orthogonal polynomial for the potential V. Pulling out the 
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integrals of the determinant leads us to 

PnNiZo) = {-IT 7^^ 

(n;=i^-') det (.^-^)"^^^e-^i:L.^(^.)d2". 
^ lur. E.e5. sign(a) (n;^. ^-^Q.;) det {z})l^,e--^^n-^ d^^z 

Inn E.e5. «ign(^) (n;.i ^to-;) det e-^E?.. v(.,) d2n, 

n 

1.3. The potential V{z) = In this case, we can remove the cut-off by setting D = 

C and immediately find that the orthogonal polynomials are nothing but the monomials. 
So we get exphcitly 

N-l 



^ n=0 'J 

Considering now the limit N ^ oo, we get for the eigenvalue density 



lim = — lim I x^-^e-^'da; 

if 

*0 



At'-^to a/"" ''^ ^"^ TT^o ./l£li 



— lim r le-^(--i-^«s-)da; 

Trio ■'v-^-oo V 2tt jMi x 



*0 



1 /iV /-^ 1 



(10) 



Trio ^^0° V 27r 1 + x 1 ^7*0' 

In the continuum limit N ^ 00 therefore, the eigenvalues uniformly fill the disc with 
radius -\/io- 

Proposition 1.8. Choosing the radius A = we find 





\z\ 




\z\ 


\ 27rto' 




\z\ 


I 



(11) A,(n,A,o)=fns.(x)) n^En#' 



w/iere Ej(a;) = e *o 

Proof. Since the orthonormal polynomials (gn,Af)$^Lo) <ln,N{z) = ./ '^S+i z^, for the po- 
tential on C are still orthogonal on .8^(0), namely 

Qi,Niz)qj,N{z)e *oi d z^Sij- / qi,N{z)qj,N{z)e *oi d z^SijEi+i{x), 

C\Ba(0) Bx(0) 
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we get 



n\{N 



1 /"e-ej:£.l^.l>etta_,fe)).':;.,rd-. 



Ba(0)"x{C\Ba(0))^ 



1 

^ ' o-eSjv «=1 j=n+l 

.j=l / n<...<i„j=l ^ 

Proposition 1.9. T/ie probability of finding exactly n eigenvalues in a disc of radius 
A = around any point a in the interior of the disc is asymptotically for N —>■ oo 
given by A^in, X, 0) . 



Proof. According to Proposition 11.41 the probability A^in, A, a) is given by 

N-n , .J. 

where the correlation functions R^^\ ^ ^ k < N , can be determined by 
with 



Kn{w,z) = NkN{wz)e ^^o^' ' \ kN^z) = — 2_^— re 'o . 



m=0 

Now kN{\z\^) = ^R^^\z) for all z G C, and so calculation (ITOi) tells us that, for z in some 
domain containing the interval (0,to); we locally uniformly have k]^[z) = ^ + o(e~^^) 
in the limit N ^ oo for some e > 0. Therefore, there exists some e > such that in the 
limit N ^ oo locally uniformly in w and z 

K]^[a + Xw,a + Xz) = N \ h o(e ""j e ^'o^' ' ' ' *" 



TXt 







So, for 1 < A; < AT, 



-^i?!;)((a + A.OtO = (7-^ + o(e--f)) j:sign(a)e-t^?-^-^(— «) 

is asymptotically independent of a and so is the probability AN{n, A, a). □ 

2. Polynomial curves 

2.1. Harmonic moments. 

Definition 2.1. A polynomial curve of degree d is a smooth simple closed curve in the 
complex plane with a parametrisation h : G 'C ^ <C of the form 

d 

(12) h{w) = rw + ''^^ajW~\ \w\ = 1, 

j=0 



RANDOM NORMAL MATRICES AND POLYNOMIAL CURVES 



9 



with r > and 7^ 0. The standard (counterclockwise) orientation of the circle induces 
an orientation on the curve. We say that a polynomial curve is positively oriented if this 
orientation is counterclockwise, that is if the tangent vector to the curve makes one full 
turn in counterclockwise direction as we go around the unit circle. 

Proposition 2.2. Let •y be a positively oriented polynomial curve with parametrisation h 
of the form ffT^ . Then h, viewed as homolorphic map on C^, restricts to a biholomorphic 
map from the exterior of the unit disc onto the exterior of y. 

Proof. We have to show that h'{w) 7^ for all w in the complement of the unit disc. Let 
if: denote the tangent vector map w ^ t{w) = h'{w)iw = i{rw — ^ jajW"-''). Since 7 is a 
simple closed curve, the map w h-^ t{w)/\t{w) \ is a map of degree one from the unit circle 
to itself. Therefore, we have 

If If t'{w) , 1 / t'iw) 



— ^ dargftfw)) = — ffi — — -A.w = N^ <h — — - dw. 

Here > denotes the number of zeros of t[w\ counted with multiplicity, in the com- 
plement of the unit disc and R is so large that it contains them all. The latter integral 
is one as can be seen by sending R to infinity. Thus, = 0, and h! has no zeros in the 
complement of the unit disc. □ 

A simple consequence of this proposition is that a polynomial curve is uniquely para- 
metrised by a map of the form ffT^ with r > 0. Indeed, any other conformal mapping 
of the complement differs by an automorphism of the complement of the unit disc. But 
non-trivial automorphisms are given by fractional linear transformations which do not 
preserve the conditions. 

In the following, we therefore mean by the parametrisation of a polynomial curve always 
the parametrisation of the form ( |T2l) . 

Definition 2.3. The exterior harmonic moments {tk)kLi of a polynomial curve 7 encir- 
cling the origin are defined by 

tk = (p zz~^ dz, e N. 



27riA; 

The interior harmonic moments {vk)'^^i of 'j are 

Vk = — (f zz^ dz, k e¥\. 
27ri 

Proposition 2.4. Let 'j be a positively oriented polynomial curve of degree d encircling 
the origin. 

(i) The exterior harmonic moments tk of 'j vanish for all k > d + 1. 

(ii) There exist universal polynomials Pj^k G 2;[r, Oq, . . . , aj^k], ^ ^ k < j , so that for 
l<k<d+l, 

d 

(13) ktk = dk-ir^''^^ + ^ ajr^^Pj^kir, oq, . . . , aj_k)- 

j=k 

Moreover, Pj k is a homogeneous polynomial of degree j — k + 1 and it is also 
weighted homogeneous of degree j — k + 1 for the assignment deg(aj) = j + 1, 
deg(r) = 0. 
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(iii) The area of the domain enclosed by 7 is nto where 

d 

(14) t,=r'-Y,j\<^j?- 

Proof. Let h be the parametrisation of the polynomial curve 7. Then, since 7 encircles the 
origin, h{w) 7^ for |w| > 1. Hence, the contour in the formula for may be computed 
by taking residues at infinity. For > 1, 

ktk = — (p h{w^^)h' {w)h{wy^ dw 
2vri 



The integrals in this sum vanish if j < k — 2. The formula for tk in terms of r and (oj)^^^ 
is obtained by expanding the geometric series and picking the coefficient of w~^ in the 
integrand. This proves (i) and the first part of (ii). The homogeneity property is clear. 
The weighted homogeneity follows by rescaling w in the integral. The same formula can 
be used to compute to, but the first term rw~^ in h{w~^), which does not contribute to 
the integral and was omitted for k > 1, must be added here. □ 

Theorem 2.5. Let t2, . . . ,td+i be complex numbers so that \t2\ < |. Then there exists a 
6 > so that for all to, ti with < to < 6 and < to(| ~ l^2|); there exists a unique 
positively oriented polynomial curve of degree less than or equal to d encircling the origin 
with area nto and exterior harmonic moments (tfc)^i with tk = for k > d + 1. 

Proof. The idea is to invert the map (r, oq, . . . , ad) ^— > (to, • • • , td+i) defined by ( fT3l) and 
ffT4l) for small r and oq. Set aj = r~^aj, p = r"^ and consider the resulting polynomial 
map 

F : (p, ao, • • • , Oid) (to, • • • , Wi), 
as a map from M x C""*"^ to itself. The first claim is that this map has a smooth inverse in 
some neighbourhood of any point (tg, • • • , td+i) G R x C"'"''^ with tg = ti = and \t2\ 7^ ^. 
By Proposition 12. 4^ this map is given by 

d 



ktk = + ^ ajPj^kir, ao, rai, 



j=k 



= afc_i + ^ ajPj^kip, "0, «i, • • • , ttj-fc), 1 < k < d+ 1 

j=k 

where the universal polynomials Pj^k, ^ ^ k < j , are given by 
-Pj,fc(p, «o, • • • = 



1 
27ri 



-k 
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for any sufficiently large R. By computing the residue at infinity, we can calculate Pj^k 
and thus tk up to terms of at least second order in ao and p, 

to = (1 - laiHp H , 

ktk = OLk-i - kakOo - {k + l)ak+iaip H , 1 < k < d+ 1. 

Hence, F(0, 0, 2 ^2, • • • , (c^ + 1) id+i) = (0, 0,t2, ■ ■ ■ , td+i) and the tangent map at this point 
sends (p, do, • • • , cxd) to (to, ... , id+i) with 

io = (l-4|t2np, 

kik = dfc-i - k{k + l)tfc+ido - 2{k + 1)(A; + 2)4+2^2^ l<k<d+l. 

The tangent map is invertible if \t2\ 7^ |. By the inverse function theorem, F has a 
smooth inverse on some neighbourhood of (0,0,^2, . . . ,td+i)- If 1^21 < \, F preserves the 
positivity of the ffist coordinate. 

In terms of the original variables, this means that given any (to, . . . ,td+i) with small 
to > 0, ti and such that |t2| 7^ |, there is a curve w ^— > h{w) with 

d 

h{w) = rw + r^ajw^^ and aj = (j + 1) t,+i + (9(r^). 
i=o 

There remains to show that if r > is small enough, h parametrises a positively oriented 
simple closed curve containing the origin. We ffist show that h is an immersion. Since 
h'{w) = r — raiw""^ + C(r^) and limr^oai = 2t2, we see that as long as |t2| 7^ \, h'{w) 
does not vanish on the unit circle. Similarly, we show that /i : S*^ — > C is injective: we 
have 

\h{w) - h{w')\^ = r\w - w' + 2i2{w-^ - w'~^)\ + 0{r^) 
= r\w — w' + 2t2(w — w)\ + (9(r^). 

But the expression in the absolute value can only vanish for w ^ w' if \2t2\ = 1, which is 
excluded by the hypothesis. Moreover, h{w) = rw + ti + 2rt2W~"'^ + 0{r^). Therefore, h 
parametrises a perturbation of an ellipse centered at ti. The condition on ti is a sufficient 
condition for this ellipse to contain the origin. □ 

Definition 2.6. We denote the space of all harmonic moments t = (t,t)^Q (including the 
area to of the interior domain in units of n) fulfilling the conditions of Theorem \2.5\ for 
deN by Td. 

In the following, we will always consider a polynomial curve to be positively oriented. 
2.2. The Schwarz function. 

Definition 2.7. The Schwarz function of a polynomial curve 7 is defined as the analytic 
continuation (in a neighbourhood of the curve) of the function S{z) = z on '-f. The 
Schwarz reflection p for the polynomial curve in this domain is the anti-holomorphic map 
p{z) = S{z). 

Definition 2.8. Let 'y be a polynomial curve and h its parametrisation. The critical 
radius R of ■y is then defined as 

R = max{\w\ \ h'{w) = 0, w G C}, 

which by Proposition \2.2\ is less than one. 
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The following proposition now shows that the Schwarz function is indeed well defined. 

Proposition 2.9. Let 'j be a polynomial curve parametrised by h and R its critical ra- 
dius. Then the Schwarz function S and the Schwarz reflection p of the curve 7 restricted 
to h{Bji-i \ Bpij, where Br denotes the open disc with radius R around zero, are biholo- 
morphic and anti-biholomorphic maps respectively. They are given by 



(15) S{z) = h , ^. . and p{z) = h 



1 



h-\z) J ' \h-\z) 

where h"^ denotes the inverse of h\c\Bii- 

Therefore, p maps 7 identically on itself, h{Bi \ Bji) to h{Bji-i \ Bi) and vice versa. 
Additionally, we have that p o p is the identity. 

Proof. By definition of the Schwarz function, in a neighbourhood of |m7| = 1, 

S{h{w)) = h{w~^). 
Because for w E (C\ B^ the function h is biholomorphic, we may write 



^W = ^ rrTTT > p{z) = S{z) = h 



h-\z) y ' ^ ^ ^ \h-\z 

Taking the derivatives, we find that they do not vanish for z G h{Bji-i \ Br), and so S 
and p are biholomorphic functions on H^Br-i \ Br). □ 

Proposition 2.10. The Schwarz function of a polynomial curve 7 of degree d encircling 
the origin, with parametrisation h and critical radius R, is analytic on h{<C\BR) and has 
the Laurent representation 

ci+l 00 

(16) S{z) = ktkz'-' + j + Y. 

k=l k=l 

around infinity, where (tfc)fcLi are the exterior harmonic moments of 7, {vk)'^=i are the 
interior harmonic moments of'y, and ttIq denotes the area of the interior domain of'y. 

Proof. With Cauchy's integral formula, this follows immediately from the definition of 
the harmonic moments: 

Vk = — I z^Siz) Az for G N. □ 
2vri 

Proposition 2.11. Inside a polynomial curve 7 the Schwarz function S of ^ becomes a 
multivalued function which has no singularities, unless the curve is a circle where we have 
S{z) = ^. 

Proof. Because of the representation (fTSl) for the Schwarz function, a branch of the 
Schwarz function may only diverge in z G C if z G {h{{]) , h{oo)} . If the polynomial 
curve is not a circle, we have /i(0) = 00 = h{oo). So, the only singularity is located at 
infinity. □ 
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2.3. The Toda lattice hierarchy. 

Definition 2.12. Let G N. The Toda lattice hierarchy is defined by the Lax-Sato 
equations 





N 


dtk 


~ h 




N 







(17) ^ = f[Mi?\^.], If = ^[<'^-- 

dtk ^0 



18) ^ = f[L.,M(f^], ^ = ^[L.,M(f)], 



G N, for the difference operators 



LN{to,t,t) = rN(to,t,t)eN'''o + ^aj^N(to,t,t)e~ ^ ""'o and 

j=0 



LN{to,t,t) =fN{to,t,t)e ^^'0 + ^aj^N{to,t,t)e 

j=0 



oo 



t = (tfc)fcLi; i = (^fc)fcLi; cn the space of all analytic functions in to, where 

(19) = {L%)++^-{L%)o and M^f^ = + i(L^)o, A: G N. 

i/ere 0+, Oq, and 0„ denote the positive, constant and negative part of the operator O 

tQ 

in the shift operator e^v *» . 

Proposition 2.13. Let L^q and L^q he a solution of the Toda lattice hierarchy. Then, 
with Mjf^ and M^^^ given by (|T^ . the compatibility relations 

(20) ^ - ^ = f |mW, m'^>], §^ - ^ = ^lMi?M'], 







N 


dtk 


dt, 


~ ~Q 




dMi^^ 


N 


dik 


dte 


~ To 



(21) ^ + ^ = ^-[m^^\m^;:\ 

fc, £ G N, for the equations ffTTj) anc? ffT^ are fulfilled. 

Proof. Let A;,£ G N, and let us introduce M"^^^^ = L% - M^^^ and M^J-'"'^ = L% - M^^l 
Then 

dlir dL% AT,, N 



Otk Ot£ to to 



= ^[Mi;\Mii^] + ^[Mf\Ml}\ 
to ^0 

Adding now the positive and half of the zeroth part of this equation gives us the first of 
the relations fl20l) . On the same way, we get by taking the negative and zeroth part of 



dL% dL% _ AT ^^(,) ^ Mcii)^Mf^ ^ N^Mik) ^ ^cik)^ ^cie)^ 
N^,~jn ,~jk^. iV.,~ 



dtk dti to to 

^0 ^0 



14 PETER ELBAU 

the second relation of fl20l). Now 



a* a, \ att at, «o' " " 

and therefore, since 







dtk dte to' ' V„ 2 V 54 /o 2 \ dt 

dM^S^ dMll^^ N^^.,, , „ 
^ TT^ + — M^\M' 



equation ( l2Ti) is fulfilled, too. □ 

Taking formally the limit — > cxo of the Toda lattice hierarchy, the shift operator 
eiv^*o will be replaced by a variable w and the scaled commutator ■] becomes a Poisson 
bracket with respect to the canonical variables log if and to- This leads us to the following 
definition of the dispersionless Toda lattice hierarchy. 

Definition 2.14. The dispersionless Toda lattice hierarchy is given by the equation system 
Oz Oz 

dz r dz , ~ , 

■^ = {z,Mk}, -^ = {z,Mk}, 

otk dtk 

G N, for functions z and z of w, to, t = (tk)'^=i, and t = (tfc)fcLi of the form 

oo 

z(w,to,t,t) = r{to,t,t)w + ''^^aj{to,t,t)w~^ and 

j=0 



oo 



z{w,to,t,t) = r{to,t,t)w ^ + y^aj(to,t, t)w^ . 

Here, with /+, fo, and /_ denoting the positive, constant and negative part of a function 
f considered as power series in w, 

(22) Mk = {z'')+ + ^{z% and = + l(5'=)o, G N, 

and the Poisson bracket is defined as 

ff g\ = w—— - w— — 

dw dto dto dw 

The compatibility relations for the dispersionless Toda lattice hierarchy are also given 
as the hmit N ^ oo of the compatibility relations (^01) and (1^ of the dispersionful Toda 
lattice hierarchy. 
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Proposition 2.15. Let z and z be a solution of the dispersionless Toda lattice hierarchy. 
Then and given by fulfil the equations 

dMe dMk r.. ... dMe dMk , ~ ~ , 

___ = {M„M.}, ^-^ = {M„A4}, 

dtk dti 

Proof. The proof follows exactly the lines of the proof of Proposition 12. 131 We only have to 
adapt the notation and substitute the scaled commutators ^f-, ■] by the Poisson brackets 

{-,■}• " □ 

2.4. Integrable hierarchy of conformal maps. We consider the parametrisation h of 
a polynomial curve as function of the parameter w G C^, of the area to (in units of tt) 
and the harmonic moments (tfc)fcli of the curve. Set z{w,t) = h{w) and z{w,t) = h{w~^), 
where t = (tfc)^o ^ ^ some d G N and w G . 

Proposition 2.16. On x Td, we have {z,z} = 1, where the Poisson bracket {■,■} is 
defined as 

df , dg df ^dg 

(23) {f,9}{w,t) = w—{w,t) — {w,t) -w—{w,t) — {w,t). 

Proof. Let S{-,t) be the Schwarz function of the polynomial curve defined by the set 
t G of harmonic moments. Then, for \w\ > 1, z{w,t) = S{z{w,t),t), and we get 

dz f dS dS dz 

{z,z}{w,t) = w—{w,t) Tr-iz{w,t),t) + -T^iz{w,t),t) — {w,t) 
aw \oto oz oto 

dz .dS., . .dz . , 

whose Laurent series in w around infinity has the form 1 + 0{w~^)^ as can be seen from 
Proposition 12.101 

On the other hand, using 2:(w,t) = 5'(z(w,t),t) for |w| > 1, where we defined 5(2;, t) = 
^(z, t), we get 

dz dS 
{z,z}{w,t) = -w—{w,t) — {z{w,t),t), 

which has a Laurent series of the form 1 + 0{w). 

So, by analytic continuation in w, we have {z, 5} = 1 on x 7^. □ 

Proposition 2.17. There exists a function Q : {{z,t) G C x 7^ | 2 G (Dj)_} C with 

dfl dQ 

(24) S{z{w,t),t) = —-{z{w,t),t) and \ogw = ——{z{w,t),t) for \w\ > 1, 

dz dto 

where S{-,t) denotes the Schwarz function of the curve 7* defined by the harmonic mo- 
ments t&Td and {Dt}^ is the exterior domain of the curve 7f . 

Additionally, every such function has an asymptotic expansion of the form 

d+l 00 , X 

(25) n{z, t) = Y, tkz' + to log z - -vo{t) - J2 

k=l k=l 
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around z = oo, where (ffc(t))^i are the interior harmonic moments of -yt and, with 
z{w, t) = r{t)w + 0{1) for w oo, 

(26) |^(t)=21og(ra)). 

Proof. For \w\ > 1 and t & Td, we write z{w,t) = S{z{w,t),t) and get with Proposition 
12.161 the compatibihty relation 

dS , , , , dz , , dS . , , , dz , , 
Q^{z{n^.tXt:) = ^^{w,t) - -{ziw,t),t) — {w,t) 



^( A Bi^^t) dz 1 



for the equation system ( 12^ . which imphes the existence of such a function Q. 

For the asymptotic expansion of Q, we integrate the expansion ( fT6l) of the Schwarz 
function S with respect to z and get with the integration constant — |fo(^) the equation 
(l25l) . Equation (1261) is now imphed by the zeroth order in w of the relation 



logw; = g^^{z{w,t),t) = logw + log(r(i)) - 2^(^) + (^(^"')- ^ 

Proposition 2.18. On x %, we have for 1 < k < d + 1 

dz dz dz ~ dz ~ 

where the Poisson bracket {■,■} 2S defined by ( l23l) anc?, wi/i i/ie notation /+, /_ anc? /o 
/or t/ie positive, the negative, and the constant part of the Laurent series of f in w, 

(28) Mk{w,t) = {z'{w,t))+ + ^{z\w,t))o, Mk{w,t) = {z\w,t))^ + ^{z'{w,t))o. 
Proof. We choose a function Q as in Proposition 12.171 and define for \w\ > 1 and t&Td 

Mk{w,t) = ^{z{w,t),t). 

Considering now only \w\ > 1, we get, using ^{z{w,t),t) = logw, 

dz d^Q 

dz f d f dVt \ d'^Q dz 

dz , . 
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Similarly, we have, using additionally z{w,f) = ^{z{w,t),t) and {z,z} = 1, 

d'^il ^ dz d'^Vl 

{Mfc,z}(w;,t) = ■g^g^{z{w,t),t) {z,z}{w,t) - w—{w,t)-^^-^{z{w,t),t) 

d f dVt \ dz d'^Vt 

dz ( d ( dVt \ dz d'^Q 

dz , . 

Solving these two equations for we find, again using {z, z) = 1, 

dMj, dz dz dz dz 



dw dw dtk dtk dw 



With z{w,t) = S{z{w,t),t), where as before S{z,t) = S{z,t), we then have 
'dMk, A f dz , ..dS 



Therefore, substituting the asymptotic expansion flzoj) for Vt 



M,{w,t) = {z\w,t))^ + {z\w,t))o - l^^it). 

Considering now the derivative by to, we find as before 

dMk , . dz , ^dz , , dz , . dz . . 

Writing = rii^w^"^ + 0{1) for if ^ 0, we see 

dMk, - 1 dr , . , 1 (9^fo , . X 

--^ t = — a + 0{w) = ^j^{t) + 0{w). 
dto r\t) dtk 2 dtodtk 

Therefore, using the freedom of choosing the to-independent part of Vq, we achieve 

Mk{w,t) = {z\w,t))+ + ^{z\w,t))o. 

By analytic continuation in w, we see that Mk, defined by ( l28i) . fulfils the relations ( 1271) 
for all w e C and teTa. 
Setting now 

Mk{w,t)=Mk{w-\t), 
and using z{w,t) = z{w~^,t), we get the corresponding relations for Mk- □ 

Therefore, the functions z and solution of the dispersionless Toda lattice hierar- 

chy in the first d+1 time variables {tk)k=i fulfilling additionally the so called dispersionless 
string equation {z, z} = 1. 

Furthermore, because of the additional restrictions tk = tk, I < k < d + 1, and 
z{w,t) = z{w~^,t), {w,t) G X Td, compared to Definition 12.141 of the dispersionless 
Toda lattice hierarchy, the equations for the time variables {ik)k=i are equivalent to those 
for the variables {tk)k=\- 
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3. The eigenvalue density in the continuum limit 

3.1. The equilibrium measure. We now turn our attention to the hmit iV — > cxo. To 
describe a configuration in this hmit, it is useful to pass from the sequence of eigenvalues 
to measures. So let us introduce the point measures 

1 ^ 

Sz{A) = Xa{z), zeC,AGC, and S.^-J^S^^z^ (^^)^=l C C, 

i=l 

where xa denotes the characteristic function of the set A G C Then, in the domain Dq , 
the probability distribution Pjv may be written as 

where we introduced the functional 

= / Viz) dfliz) + [[ log\z- C\-' d/i(C) dfxiz) 
J J Jz^C 

on the space of all Borel probability measures. (Again, we will drop the index indicating 
which potential we use.) 

Lemma 3.1. The space of all Borel probability measures on D is sequentially compact. 

Proof. By the theorem of Riesz-Markov each Borel measure jionD corresponds to exactly 
one positive, linear functional (j)^ e C{D)* and by the theorem of Alaoglu, the closed 
unit-sphere in C{D)* is weak-*-compact. Therefore, for each sequence (yU„)^^ of Borel 
probabifity measures, the sequence {(t>iJ.n)'n=i contains a weak-*- convergent subsequence 
(</'Mn(fc))r=i> that is 

30gC(Z})* : 0,^^^,(/)^0(/) V/eC(D). 
Now we find a measure on D with = 0^. This measure fulfils 

J /d//„(fc) ^ J fdfi (k ^ oo) V/ e CiD), 

and hence, is again a Borel probability measure. □ 

Lemma 3.2. Let fi and jl be Borel probability measures such that the function log — (^|~^ 
is integrable with respect to ii® ii and to ji® ji. Then log 1^; — is also integrable with 
respect to ^® fi. 

Additionally, we have the inequality 



(29) j J log \z - C\-' d(/i - /.)(C) d{jl - ii){z) > 

with equality if and only if /i — jl. 

Proof. We start with the distributional identity 

/ log \z\-'^Aip{z) d^z = -2tt^{0) 

for any smooth function (/? : C — > R with compact support. Introducing the Fourier 
transform 



^(k) = — [ (^(^)ei('=^+^^") d'^z, keC, 
27r 
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of 93, we find 



/ \og\z\-^/\v{z)d?z = - [ -^{\k\'0{k))d'k 



|2 

where / : C — >■ [0, 1] denotes a continuous function which is one in the vicinity of zero 
and becomes zero at infinity. So, we have for all z G the equation 

= i X ^ ("^^''^''^ - ^(^)) + ^(^) 

for some real constant C{f). 

Therefore, with Tonelli's theorem, we see that 

\og\z-C\-U{~^-^^md{~^-^^){z) = ^ |e^('=^+^^")d(/2-/i)(z)'d2A; 

is non-negative and finite, which immediately implies the integrability of log|2; — 
with respect to jj,® jl. To achieve equality in (1^ . we need 

^l{kz+-k-z) ^^^^^ ^ j ^^(kz+k-z) ^^^^^ 

for all A; G C, which reads fi = fi- □ 

Definition 3.3. A Borel probability measure Hq on D without point masses is called an 
equilibrium measure for the potential V on D G C if 

/(/io) = inf /(/x), 

tJ.&M{D) 

where A4{D) denotes the set of all Borel probability measures on D without point masses. 
We then set Iq = I{fio). 

Since I{j^^^Xd) < 00, where denotes the Lebesgue measure on D, Iq is finite. 

Theorem 3.4. Every continuous potential V on a compact domain D has a unique equi- 
librium measure. 

Proof. To show the infimum is achieved, we choose a sequence in M.{D) with 

li.lJ'n) h- Because of Lemma ISTTl there exists a convergent subsequence {nn{k))^=i of 
(/in)$^i and a Borel probability measure /i with /i„(fc) — > /z. 

To prove that /(/i) = Jq, we estimate with an arbitrary real constant L 

lim /(/i„(fc)) = lim / V {z) d^n(k){z) + lim // \og\z - C\~'' d^n{k){() '^^^n(k){z) 

fc— >oo fc— *oo J k^oo J J 

> j V{z) dfi{z) + lirn^ jj min{log|2;-Cr\^}d/i„(fe)(C)d;U„(fc)(z). 

Approximating uniformly the second integrand according to the theorem of Stone- Weier- 
straB up to some e > with a polynomial in z, z, ( and ( and using Fubini's theorem, 
we get 

/o = J-im > j V{z) dfi{z) + jj min{log 1^; - ^} d;u(C) dn{z)-2e. 



27r 
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Letting first e and then L —>■ oo shows that fi has no point masses (otherwise the 
right hand side would diverge) and I{fi) = Iq. 

Next we want to show that there is exactly one measure fi G Ai{D) with /(/i) = Iq. So 
suppose fl G Ai{D) fulfils I (ft) = Iq, too. Then we consider the family 

Ht = tjl + {I — t)fi = /i + t{jl — /i), t G [0, 1], 

in M{D), expand with regard to Lemma [3.21 the functional lifJ^t), and obtain 



/(/i,) = /(/x) + t / (V{z) + 2 / \og\z-C\-'df,{0] d{j2-fi){z) 
(30) J \ J J 

+ t' [[ log\z-C\-'d{fx-f^){C)difi-fi){z). 



Lemma 13.21 now states that the coefficient of t"^ is non-negative, and so the function 
t 1—^ I{fJ't) is convex on [0, 1]. In particular, for all t G [0, 1], 

I{fit)<tI{jl) + {l-t)I{fi) = Io, 

which implies I{fit) = for all t G [0,1]. This requires the last summand in (IHUj) to 
vanish, and so, again with Lemma [3. 2 [ we see that /i = /i. □ 

Proposition 3.5. The probability measure fi is the equilibrium measure for the potential 
V on the domain D if and only if the function 



(31) E{z) = V{z) + 2j\og\z-C\-' dfiiO 
fulfils the relation 

(32) j E{z)dji{z)> j E{z)dfi{z) for all fi e M{D) . 

For the equilibrium measure fi and Eq = J E{z) dfi{z), we additionally have 

(33) E{z) = Eq fi-almost everywhere. 

Proof. Let us first assume fi is the equilibrium measure. Then, for an arbitrary measure 
fl G Ai{D), the condition 

^/(/iOL=o^O' f^t = tfi + {l-t)fi, tG[0,l], 
has to hold. This means 

j (viz) + 2 J log\z-C\-'dfi{C?j d(/i-/i)(^)>0. 



which immediately implies equation 

Assume on the other side that fi fulfils the condition ( !32l) . Then we obtain for the 
equihbrium measure fio 

I{fio) = I{fi) + J (^V{z)+2 J log\z-C\-'MC?j d(/io-/i)(^) 

> 

and so fi = fiQ. 
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To prove relation (l33l) . we consider for the equilibrium measure the set 

B = {zeD \ E{z) < Eo}. 

If fi{B) > 0, the variational principle for the measure fl given by d/i 
yield 

Eo<l E{z)^^df,{z)<Eo, 

and so ^{B) has to vanish. Thus, we get relation 

Corollary 3.6. If for a measure fi G Ai{D) the function E defined by equation flHTl) fulfils, 
for some real constant Eq, E{z) = Eq on the support of fi and E{z) > Eq everywhere, 
then fj, is the equilibrium measure. 

We may give this variational principle an electrostatic meaning. We consider the equi- 
librium measure fi as density of charged particles in the plane with total charge 1. Then, 
since 21og|z|~^ is the Green's function in two dimension, the function E becomes the 
electrostatic potential produced by the external potential V and the charged particles. 
Corollary 13.61 now tells us that the electrostatic potential is constant inside the charged 
area and increases outwards, so the configuration described by the measure fi is stable. 
The functional I becomes in this setting the electrostatic potential energy of the config- 
uration, which by construction is minimised by fi. 

Out of this analogy, we see that for an absolutely continuous measure fi with dfi{z) = 
p{z) d'^z, the density p will solve the Poisson equation AV^ = —iiip inside the support of 
p. 

3.2. The eigenvalue density. Let us now calculate the eigenvalue density and, more 

(k) 

generally, the correlation functions R)^ , ^ < k < N, out of the equilibrium measure p. 
To this end, we first consider the behaviour of the probability distribution around the 
equilibrium measure. So consider the sets 

AN,r, = {ze I < Jo + r]}, 7f> 0. 

Let us further assume that the boundary dD of the cut-off D is a twice continuously 
different iable curve. 

Lemma 3.7. For alle > 0, there exists an absolutely continuous Borel probability measure 
Pe on D so that d/ie(z) = Pe{,z) d?z with p^ G C{D) and 

Hl^e) <Io + e. 

Proof. Set B^{z) = B^{z) fl D. Since the boundary of D is chosen as twice continuously 
different iable curve, we find an eo > such that the circle dBs{z) cuts the boundary dD 
of D for all z E D and e G (0,eo) at most twice. Then, for any 6 G (0, |), we find an 
El G (OjEq) such that the area Kei^) = \{B^{z)) of B^{z) fulfils 

5 < < 1 for aWzeD and e G (0, ei). 



Xb 



dp would 



□ 
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Let now fi be the equilibrium measure and set Pe{z) = Then, for all 

continuous functions ip G C{D), 



D 



D 



< 



< 



D 



max 



d/i(C) 



So the measure /i^, defined by dfj,e{z) = Pe{z) d'^z, converges weakly to the equilibrium 
measure p. Moreover, we have 



log|z — w| ^p^{z)pe{w) d'^zd'^w — / / logiC — ^ dp{() dp{uj 



D2 



< 



BP (Ox BP (u) fie[Z)Ke[W) 



log 





-1 


z — w 


d^zd^w 







d'^zd'^w-1 



' BP (Ox BP (c.) l^e{z)Ke{w) 

Now the second term tends to zero in the limit e ^ 0, and because of 

-1 

d^^d^u; 



dp{C) dpiuj) 
|log|C-^ll dp{C)dp{uj). 



1 



BP(OxBP(uj) t^e\Z)l^e\W) 



log 



z — w 



< 



< 



I log 11 + 2; — It'll d'^zd'^w 



|log|l + t;|| d\! < C 



log 



1 + 



2e 



IC-^I 



(£^0), 



for some constant C > 0, the first term, too, vanishes in the limit e ^ 0. Therefore, 
I{Pe) ~^ h for s ^ 0. Hence, we find for any e > some e > such that the absolutely 
continuous measure p^ = -pjjxjPe fulfils /(/Xe) < /o + □ 

Lemma 3.8. The probability P]\f(D^\AN^^) drops for N oo to zero faster than e"^^^''. 

Proof. Using Lemma [3.71 we choose a Borel probability measure /i^ with 

Hl^v) < Io + ^ and dpr,{z) = pr,{z) d^z, pr, E C{D). 

Then, with Jensen's theorem and / /(^z) H dpr^^Zi) < /q + | + o(l) for — oo, we get 



Zn> 



N 



g-7V2/{5,)-ElLilogp^(2,) 



i=l 



and therefore, 

JL 



N 



,N-^(lQ+l)+o(N^)-N^(lQ+r,) 



^0 



Zi = o e 2 



□ 



i=l 



Lemma 3.9. For rj > 0, any sequence {zj^)n£¥i, z^ E A^^ri, and any convergent subse- 
quence (52^j^j)neiN of the corresponding 6 -measures, we have 

h < H^v) < h + V, where i/^ = lim 6, e M{D). 
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Proof. Using Zjv(n) ^ ^N{n),ri, we obtain with the cut-off L e R: 

Io + V> Yl + Yl min{log |z7V(n),i - ZNin),j\~\ L} 

= I V{z)d5,^^Jz) + llrmn{log\z-C\-\L}d5,^^jC)d5,^^Jz)- 



N{n) 



Sending first n and then L to infinity brings us to G Ai{D), and therefore Iq < li^q) < 
Jo + r^. □ 

Theorem 3.10. For all if G C{D^), we have the equality 

k k 

(34) jhn [ ^^((^.)tO<'^((^^)ti) n d'^^ = / n 

^'^■^D'^ i=l J i=i 

That is the measure T^-Ryv'* ((-2i)i=i) Yli=i ^'^^i on converges weakly to nf=i 

Proof. Substituting in the left-hand-side of equation the definition of the correlation 
functions and turning our attention on the highest order in A^, we obtain (because P/v is 
invariant under the symmetric group) 

^^^kR'^) = ^kj^^^i^)R'^i^)^''-- 

N 

= ]^. E j^^v{{z^,)U)PN{z)d?'^z + o{l). 



N 



Because of Lemma 13^ we may now integrate over the set AN,r) instead of D . So let 
the continuous function ■^J2f{i^i])j=i) ^^^^ maximum on the compact set A^^^ at 
(, and set u^^r] = ^c- Then 

«lv.«fe = l « = 1 

Because of Lemma [XTl we find a convergent subsequence i'N{n),r} {n ^ oo) with 



Now, because of Lemma 13.91 letting ^ 0, a subsequence of Ur, converges to the 
equilibrium measure /x, and thus, 



Arguing on the same way for the limes inferior concludes the proof. □ 
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4. The equilibrium measure for a polynomial curve 
4.1. The general result. 

Theorem 4.1. Let d G N. Then, for any set C C with ti = 0, |t2| < \, and 

tk = for k > d + 1 and any compact domain D G C containing the origin as interior 
point such that the function 

d+1 



U{z) 



- 2 Re ^ tki 



k=2 



is positive on D\ {0}, there exists a 5 > Q so that for all < to < 6, we have (tfc)fcLo ^ "^d, 
and the equilibrium measure for the potential V = -^U on D is given by 

1 

where denotes the interior domain of the polynomial curve 7 defined by the harmonic 
moments {tkj'^Lo! ^'^^ is the Lebesgue measure on D^. 

Therefore, the normal matrix model solves the problem of finding the interior harmonic 
moments {vk)'k'=i out of the exterior harmonic moments {tkj'^Lo-, namely (as can be seen 
from Theorem 14.11 and I3.10p 



Vk 



lim , 



Lemma 4.2. Let 



V{z) = i 



riD) 



d+1 



ti{M^)VN{M) dM. 



\z?- 



2 Re ^ tk'. 



k=l 



and denote by 7 the polynomial curve defined by the harmonic moments (tjfc)^Q. Then, 
in the interior domain Dj^ of^, the function 

E{z) = V{z) + 4 / log 7 - 1 ' d^C 

TTto Jd+ C 

is equal to zero, and in the exterior domain C \ -D+, its gradient reads 

1 



(35) d,E{z) = -{z-p{z)), 

to 

where p is the Schwarz reflection on 7. 

Proof. To verify the first statement, we use Green's theorem and obtain 







- / log 




Jd+ 





d'c = -kr+Re— 




z 



-1 



C + 



C-z 



dC. 



Integration by parts of the second integrand yields immediately 

-1 



2 





z 


/ log 




Jd+ 





d2( = _k|2_2Re— d) log 
2m 



and expanding the logarithm around z = leads us to E{z) = for all z G -D+. 

For the proof of the second part, we write the Schwarz function S of the curve 7 
a.s S = Si + Se, where Si is analytic in and Se in the complement C \ D^. From 
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Proposition 12.101 we already know that Si{z) = Ylkt.i^^tkz''^^. For the function 5*6, we 
find with Cauchy's integral and Stokes' formula 

1 iC-5i(C)^, If 1 ^2. 



Seiz) = -— (h^^—^dC = - -d^C, ze€\D^. 

27ri C~ z Jd+z-Q 
Therefore, for all z G C \ 

d,E{z) = ^(z-^-^) = ^{z-p{z)). □ 
to V ^ to 

4.2. The Gaussian case. In this case, where the polynomial curve is an ellipse, we are 
able to calculate the equilibrium measure on C explicitly. 

Proposition 4.3. For any to > and any ^2 G C with \t2\ < \, the equilibrium measure 
fi for the potential 

to 

on C is given by 

1 , 

6162 

where denotes the interior of the ellipse 



Re{Vhz)\ ImiVhzY /l + 2|t2|, , /l-2|t2 

(36) "2 \ 72 = \t2\, Ol 



bi bi ' \/i-2|t2r' 1/1+21^2 

Proof. As polynomial curve, the ellipse ( !36l) has the parametrisation 
(37) h{w) =r{w + 2t2W-^), r = ^i-^-^. 

We check that the given measure is the equilibrium measure by verifying the conditions 
of Corollary 13.61 To this end, let us introduce for |i(7| > 1 the function 

-1 

^2 











'^to Jd+ 





Integrating equation and its complex conjugate analogon, we get for S{w) the ex- 
pression 

g(w) = - (\h{w)\^ -\h{l)\^ -2Re r h{uj~^)h\u)do}\ . 



to \ J I / 

Substituting in this expression relation (!37j) for /;,, we obtain 

to8{w) = r^{\w + 2t2W-^\^ - |1 + 2t2r) 

- 2r^ Reii2W^ + (1 - 4\t2\'^) logw + t2W~^ - (^2 + ^2)) 
= r^{\w\^ - 1 - 4|t2|' + 4|t2|Vr' + (1 - 4|t2n logd^r')) 

+ 2r^Re{2t2Ww'^ - t2iww'^)\w\''^ - t2iww~^)\w\'^) 
= r'{\w\' - 1)(1 - 4|t2| Vr') + r^{l - 4|t2n hg{\w\-^) 

- 2r^{\w\^ - 1)(1 - \w\-^) Re{t2ww-^) 

= r'i\w\' - 1)(1 - 2|t2|)(l + 2\t2\\w\-') + r\l - A\t2\') \og{\w\-^) 
+ 2r^{\w\^ - 1)(1 - |«;r^)(|t2| - Re{t2ww-^)). 
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We are now ready to start estimating S{w) for \w\ > 1 and \t2\ < |. The last bracket we 
can estimate by 

1^2! — Re(t2ww~^) > 1^2! — \t2ww~'^\ = 0. 
There therefore remains to show 

(Iwp - 1)(1 + 2\t2\\w\-^) + (1 + 21^2!') logd^r') > 0, 

which follows immediately out of the following lemma. 

Lemma 4.4. For a G [0, 1], the function 

f{x) = (x — 1)(1 + ax^^) — (1 + a) logx 

is non-negative on the interval [l,oo). 

Proof. For the function / and its derivatives /' and /", we have 

/(1) = 0, /'(1) = 0, and /"(a;) = l(^l + a-^^ >0. □ 

So we showed that E{z) > for all z G C \ D+. Because of Lemma [4.21 we also know 
that E is zero in the interior domain , and we therefore can apply Corollary 13.61 to see 
that yU is indeed the equilibrium measure for the potential ^ on C □ 

4.3. The proof of Theorem 14. IL As in the Gaussian case, we are going to show that 
the measure /i given in the theorem fulfils the conditions of Corollary 13.61 and is therefore 
the uniquely defined equilibrium measure. 

Because Theorem 14.11 is only valid for interior domains with small area, we are going 
to consider the asymptotical behaviour ~^ 0, where the harmonic moments (t^)^^ are 
kept fixed. To catch the asymptotical behaviour of the corresponding polynomial curve 
7, let us write its parametrisation h as in the proof of Theorem 12.51 

d 

h{w) = rw + r^ajw"^ . 

j=0 

Then, as we know from the proof of Theorem 12.51 we have for r — > 

to = (l-4|t2nr2 + C(r^) and a, = {j + l)tj+i + 0{r^), < j < d. 

Lemma 4.5. The critical radius R of 'y is asymptotically constant for r 0: 

R=./\^\ + 0{r). 

Proof. The roots of the function h'{w) = r — Yl'^^ijr^ajW~^~^ are in zeroth order at 
±^/al and {{n — l)-times degenerated) at zero. □ 

We consider now for z E D and w G h"^{D\ D^) the functions 



(38) E{z) = V{z) 











T^to Jd+ 





-1 

]2 



d C and 



(39) 8{w) = E{h{w)) = -(\h{w)\'' -\h{l)\^ -2Re j h{uj~^)h! {uj) duj ] . 



to 

We already showed in Lemma [4.21 that E{z) = for all z G -D+, so the first condition 
of Corollary 13.61 is satisfied (this is essentially the way we have chosen our potential V). 

Now, again with Lemma [4.21 we see that E > in the vicinity of the curve 7, strictly 
speaking in the domain h^B^-i \ Bi), where R denotes the critical radius of 7. Indeed, 
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if we look at the connected components of the contour hnes of the function E (which are 
smooth curves in the considered domain because there dzE 7^ 0), we see that the gradient 
vector dzE always points outwards, that is into the exterior domain of the contour line. 
Therefore, the value of E on the contour lines is increasing outwards as desired. 

A bit farther from the curve, that is for < \w\ < for some e G (0, |), the function 
£ equals asymptotically the one of the corresponding ellipse h^'^\w) = rw + ao + raiw~^, 
we denote it by S^^\ Indeed, remarking that we have t^Q^ = to + 0{r^) for the area nt^Q^ 
of this ellipse and that £^^\w) = (9(r~^^), we uniformly in w obtain 

s{w) - f = -(IM^)I' - \h^°\w)\^ - + 

^0 

+ -Re/ {h{uj-')-h^''\uj-^))h'{iu)dtu 



to 

r 

to 



+ - Re / (h'iuj) - /i(°)'(lu))/i(°)(cu-^) dcu 



w) 



to 

= 0{r) + 0{r'-^') + 0{r) + 0{r''-''') 
^0 (r^O). 

Now, from the proof of Proposition 14.31 we know that there exists a C > such that 
£^^\w) > C for all w G C \ B^-i and any r > 0. Therefore, we may choose r so small 
that \£{w) - £^^\w)\ < £^^\w) and so £{w) > for all w G Br-s \ Br-i. 
There remains the domain where \w\ > . For k > 2, we obtain 



h{w)'^ - (rw)'' = ^ 



Since 



U{z) = \z\^ -t2Z^ -hz^ +oi\z\^) 

\z - 2i2z\^ + (1 - 4|t2r)kr) + od^n 



1 

2 



for 2 — i> and U > on D\ {0}, there exists a constant c > such that U{z) > c|2;p for 
all zeD. 

Therefore, for r — 0, 

V{h{w)) = V{rw) + > -Irwp + 

to 

which tends to infinity at least as r~^^. 

On the other hand, the integral over the logarithm in (1381) diverges for r ^ only as 
logr. So, for r small enough, we have E{h{w)) > for all w G h~^{D \ D^) \ B^-e. 

This proves E{z) > for all 2; G -D and therefore, with Corollary 13. 6[ that n is the 
equilibrium measure. 
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4.4. Shifting the origin. As a little generalisation, we consider the case where ti ^ 0, 
which corresponds to a shift of the origin. So we need to define the harmonic moments 
also for curves which do not encircle the origin. 

Definition 4.6. Let 

d 

-j 



h{w) = rw + ( 

he the parametrisation of a polynomial curve of degree d. Then the exterior harmonic 
moments (tfc)^^} are given by the equation system ( fT3l) . All other exterior harmonic 
moments are set to zero. 

Proposition 12.41 tells us that this definition coincides with the previous one if the origin 
is in the interior domain of the curve. 

Corollary 4.7. Let c? e N. Then, for any set (tfe)^i C C with 1^2 1 < | and tk = for 
k > d + 1 and any compact domain D G C such that the function 

d+l 

U{z) = \z\'^ -2ReJ2tkz'' 

k=l 

has exactly one absolute minimum in the interior of D, there exists a 5 > so that for 
all < to < 6 the equilibrium measure n for V = j^U on D is given by 

where denotes the interior domain of the polynomial curve 7 defined by the harmonic 
moments (tfe)fcLo- 

Proof. Let us first shift the origin by Oq, such that V{z + Oq) has its absolute minimum 
aX z = Q. Thereby, the potential gets the form 

V{z + ao) = ^ M;2p-2ReV4zM +V{ao), 
V k=2 J 

where the t'j^ are the harmonic moments of the shifted curve 7 — cto- Indeed, the coefficients 
of V and the harmonic moments depend polynomially on the shift oq. Because we know 
them to coincide as long as the origin is inside D^, they do so for all oq. 

Applying now Theorem 14. II for the shifted potential gives the desired result. □ 

5. The orthogonal polynomials corresponding to polynomial curves 

5.1. The zeros of the orthogonal polynomials. Let us consider a potential V on 
some compact domain D such that Theorem 14.11 is applicable. As in Section 13.21 we will 
assume that the boundary dD of the cut-off D is a twice continuously differentiable curve. 

Lemma 5.1. For x G (0, 1], the orthogonal polynomials {pn,N)'^=o foi" the potential V fulfil 
for all z ^ D 

lim -log|p„,jv(^)| = / log|2;-C|d/ix(C), 

where fix is the equilibrium measure for the potential . 

Here and in the following, the limit N 00, ^ —>■ x shall be understood as N,n 00 
such that 1;^ — a;| = o{j^). 



RANDOM NORMAL MATRICES AND POLYNOMIAL CURVES 29 



Proof. Using expression ([9]) for Pn,N, we get for any z ^ D 



z 

e 



ELi log{l- ^)-n / log(l- §) dM.(C)- f ) Er=i VK) ^2n^^ 



Let now w„ denote the position of the maximum of the real part of the function 

w ^ e^^=i'°sa-^)-"/i°g(i-|)dM.(C)-{^-f)ELiVK) 

on the compact set 

X X 

Then, as in the proof of Theorem I3.im 

Because of Lemma [3.91 we get in the limit N ^ oo, ^ x that 

Re (^^Pn,iv(;2)e^"^'°^(^-^)'^^-(^)^ < e°(") + o(e-^"^). 

Doing the analogue estimates also from below and for the imaginary part, we get the 
predicted equation. □ 

Lemma 5.2. In the limit N oo and — * x G (0, 1], almost all zeros of the orthogonal 
polynomials Pn.N ofV,n& Nq, lie inside the cut-off D. 

Proof. Since the function 

-log|p„,jv(^)| = / log|2;-C|d(55 

^ denoting the zeros of Pn,N, converges outside the domain D pointwise to the harmonic 

function J log I2; — CI ^l^x{,C)^ where /i^ is the equilibrium measure for the potential , 
we find that in each compact area K outside of D the ratio between zeros of pn,N lying in 
K and tends to zero. So the percentage of zeros of p„ at inside the domain D is equal 
to 

hm ^l^-^^^dz=^ls.Az)dz = l, 



N^OD,^^x 2TTm /q Pn,N\z) 2T[lXtQ 

where a is some curve encircling D and omitting all zeros of the orthogonal polynomials 
and, as in the proof of Lemma IT2l Sx,e{z) = Sx{z) — Yl'k=i ^^kZ^~^^ where denotes the 
Schwarz function of the boundary curve of the support of the equilibrium measure ix^- D 

Theorem 5.3. Let denote the zeros of the n-th orthogonal polynomial pn,N of V. 
Let further ax be the limit of a convergent subsequence of the 6-measures 6^ in the limit 

—n,N 

N —>■ 00, ^ X G (0,1]. If (Jx is supported on a tree-like graph Q, whose edges are 
smooth one- dimensional manifolds, then 

(i) the edges lie on branch cuts of the Schwarz function Sx of the curve ^x defined 

by the exterior harmonic moments (tk)'^^i and the encircled area nxto, which are 

chosen such that the discontinuity 63^ of Sx fulfils 

(40) Re{SSx{a{t))a{t)) = 0, 
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for any (local) parametrisation a of the graph Q, 

(ii) cr^ has, unless V{z) = j^lz]"^ where = 60, no point masses, and 

(iii) along the edges we have, with dz being the infinitesimal line element on the edges, 

S S ( 

(41) da^{z) = p^{z)dz, where p^{z) = ^J^^^ . 

Proof. Lemma tells us that the limiting measure ax fulfils for all z G <C\D the equation 

J log \z - CI da.(C) = 1 log \z - CI d/i.(C). 
Then, differentiating and analytic continuing this identity, gives us 

(42) I ^d(T,(C) = ^S',,e(;^), 



for all 2; G C \ ^, where, as in the proof of Lemma W?2\ Sx,e{z) = Sx{z) — Ylt=i ktkZ^~^. 

If not all exterior harmonic moments are zero, the Schwarz function Sx is according 
to Proposition 12.111 finite everywhere, and therefore the measure ax must not have any 
point masses. Setting therefore dax{z) = px{z)dz and comparing the discontinuities in 
equation (jl2D along the edges of the graph Q, we find the relation (HTI) . 

Since dax has to be a real number between zero and one, the edges of the graph for 
all segments a of the graph Q, and therefore our choice of branch cuts for the Schwarz 
function, have to fulfil the equation (HOl) . □ 

In view of the examples given in Section 15.31 it seems plausible to us to expect this 
behaviour for the orthogonal polynomials in general. 

Conjecture 5.4. The sequence {S^ t^)'n=i f^^ zeros of the n-th orthogonal polyno- 
mial pn,N of the potential V always has a convergent subsequence which tends in the limit 
N —>■ 00, ^ —>■ X & {0,1] to a measure ax which is supported on a tree like graph whose 
edges are smooth one-dimensional manifolds, except for the degenerate case V{z) = 

5.2. Connection to the Toda lattice hierarchy. We define the operators Ln, A^r and 

(k) 

Aj^ on the completion Ti. of the inner product space of all polynomials endowed with the 
scalar product (■, ■)]\f, given by equation ([8]), through their actions 

LNqn,N{z) = Zqn,N{z), ANqn,N{z) = d^qn,N{z), A^^\n,N{z) = dt^qn,N{z) 

on the orthonormal polynomials qn,N, "n. G Nq, for a potential V fulfilling the assumptions 
of Theorem 14.11 

Writing Om,n,N for the matrix elements {qm,N,0]\fqn,N)N of some operator Oj\f, we see 
that 

(43) Lm+2,n,N = 0, Am,n,N = 0, and = for m>n. 

To get the asymptotic behaviour of the orthonormal polynomials, we use Lemma [5?T] and 
will therefore again consider a cut-off D with twice continuously differentiable boundary. 

Lemma 5.5. For all m,n,N ^ Nq, we have 

\Lm,n,N\ < max|2;|. 
ze-D 

Proof. With the Cauchy-Schwarz inequality, we find 



Lm,n,N\ — \{(lm,N, Zqn,N)N\ < \/ {<lm,N , (lm,N) N \/ {zqn^N{z) , Zqn^jy) <max|2;|. □ 

V V zeD 
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Lemma 5.6. In the limit N —>■ oo, there exists an e > 0, depending only on the potential 
V and the cut-off D, such that uniformly in m,n < N 



(44) 



dD 



Proof. According to Lemma |5.1| and the proof of Theorem |4.1|, the norms {hm,N)m=o 

of 

the orthogonal polynomials {pmN)m=o ^^e potential V fulfil in the limit N ^ oo, 
f-.xG(0,l]: 



hm.M > 



D\D^ 



-mEx{z)+o{m) _ ^-mEx{0)+o{m) 



where -Dx,+ denotes the support of the equilibrium measure fix for the potential and 

Ex{z) = -V{z) +2 [\og\z-C\-' d/i,.(C). 
X J 

For the orthonormal polynomials qn.N = —7= — Pn,N, n G Nq, we therefore get in the 
limit iV 00, f ^ a; G (0, 1], f 2/ G (0, 1] 



lim — log 



,(^)g„(^)e-^^(^)d^ 



dD 



< -- mm{x{Ex{z) - i?,.(0)) + y{Eyiz) - Ey{0))), 

Z z^au 



which by the proof of Theorem 14. II is for all x, y G (0, 1] strictly negative. Moreover, since 

-1 

d^C = V{z) > 





z 


/ log 









hm x{E,{z) - E,{0)) = Viz) + ^ lim 

for all z G dD, there exists an e > such that relation (1441) holds. 
Proposition 5.7. We have the operator identity 

d+l 

(45) = -^Ai, + ktkL'},-' + Cm{D), 



□ 



k=l 



where the matrix elements Cm,n,N{D) of the operator Cn{D) tend in the limit N ^ 00 
uniformly in m,n < N to zero faster than e~^^ for some e > 0, depending only on the 
potential V and the cut-off' D. 

Proof. Writing the multiplication by z as differential operator acting on e~^^*^^-* and in- 
tegrating by parts, we get 

d+l 



{(lm,N, Zqn^N)N — -^{dz(lm,Ni (ln,N)N + ktk{z^ ^ (lm,N i Qn, 



NJN 



k=l 



1 

+ 2i 



,^(z)g„,;v(^)e-^^(^) dz, 



dD 



where the last term is by Lemma 15161 for m,n < N oi the form o(e ) for some e > 0. □ 

Let us now introduce the projection 11^ from the operators on H onto the {k-\- 1) x (A; + 1) 
matrices defined by 

nfc(O) = (Om,n)m,n=o fo^^ any operator O. 
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Corollary 5.8. The operators Ln and L*^ fulfil in the case d > 1 
(46) Il^.d{[IlN{L%),UN{Lr,)]) =^^ + 0^^,(0), 

where the matrix elements Cm,n.N{D) of the operator Cn{D) tend in the limit N —>■ oo 
uniformly in m,n < N faster to zero than e~^^ for some e > 0, depending only on the 
potential V and the cut-off D. 

Proof. Because of the relations fHHj) . we have for any operator O and any /c G N the 
identities 

n^_fe(njv(0)n^(L^)) = njv-fe(OL^) and n^(0)n^(A) = n^(OA). 

Substituting now the relation fH31) for L^, we get 

t "^^^ 

k=l 

= ^ + n^_,([n^(Cjv(/^)), n^iLr,)]). 

With CN-d{D) = Ilj^-d{[llN{CN{D)), n^(L;v)]), this gives equation 1^. □ 
Proposition 5.9. For 1 < k < d + 1, the operators Ln and Ly fulfil the equations 

^^^^ ^-^[^^ W^f^ ' 

where, with the notation 0+ for the lower triangular, 0_ for the upper triangular, and 
Oo for the diagonal part of the operator O, 

(49) M« = (L^)+ + i(L^)o, M(f)* = ((L^)'=)_ + l((L^)% 

Proof. For all n G Nq, we have 

oo r\ OO OO 

e,m=o ^ e=o ^ e,m=o 



and so %^ = [LatjA^^^]. Now, because of 



Ak) _( dqn,N\ _ [dqm,N \ N 

Am,n,N - \ qm,N, J " ^ \ ^"'^ / ~ t^^^""'^' ^ ^"'^'^ 

for all m,n & No, we have {A^^'')^ = —^{L%)_, and introducing the orthogonal polyno- 
mials Pn,N = \/hn,Nqn,N, ^ G Nq, for V, wc find 

gn,7V, = \/K^— j== = I ^1^, qn,N for all n G Nq, 

dtk J N 9tk ^/K^ \ dtk ) ^ 

which implies (^^■')o = ~^(-^Ar)o- So, by setting M^-* = L^ + ^y4^\ we get the relations 
^9]) and fulfil the first equation of (H7j) . Remarking that 

^miV = (^m,N, all m, n G No, 
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we see that also the first equation of ( H8l) and therefore, by taking the adjoint equations, 
all four equations of (H71) and (HHl) are fulfilled. □ 

Again, as in the case of the polynomial curves, only the equations for the times {tk)ft\ 
need to be considered. Moreover, equation (HUj) provides us with a dispersionful analogon 
of the string equation. 

5.3. Application. Let us for some c? G N consider the potential 

V{z) = ^{\z\' - U+i{'/+^ + ^'^+^)), Wi > 0, 
Co 

on a cut-off D chosen as disc centered in the origin such that V{z) > for all z G -D \ {0}. 
Additionally, for ci = 1, we require t2 < \. 

Using the weighted homogeneity of the universal polynomials Pj^k in equation ( |T3l) . we 
find for small to > the parametrisation /i of a polynomial curve with only non- vanishing 
harmonic moment td+i and encircled area hIq: 

(50) h{w) =rw + aw~'^, a = {d + l)td+ir'^, 

where r is the smallest positive solution (which does exist for small to) of the algebraic 
equation 

(51) to = r^-did+l)Hl^,r'''. 

Proposition 5.10. There exists ane > such that the orthonormal polynomials {qn,N)n=o 
for the potential V fulfil ( with q-n,N = for G N j for < n < N — 1 the recursion 
formula 

(52) ZQn^N = ^n+l,N(ln+l,N + Q'n-d.NQn-d.N + ^(^ '^^)(ln+l-j{d+l),N ■, 

where 

d 

(53) dn,N = {d+l)td+i\[rn+k,N + o{e-'''), < n < N - d, 

k=l 

and (setting r_n,N = for n E ¥\o) the rn,N are recursively defined by 

(54) |a„,jvr=r2+i,iv-$^|an-d+fer-^(n + l)+o(e-5^^), < n < N - d. 

k=l 

Proof. Proposition 15. 71 together with the restriction fHHj) shows us that for m,n < N only 
the matrix elements Lm,n,N with —l<n — m<d are not zero or at least exponentially 
decreasing for N oo. Moreover, the symmetry z t-^ e'^+'^z of the potential V implies 
that the orthonormal polynomials have the form 

qj(^d+i)+r,Niz) = z''qj(d+i)+r,N{z'^'"^), < r < d, j G No, 

for some polynomials qj(d+i)+r,N each with degree j. Therefore, the only non- vanishing 
matrix elements are Ln+i~j(d+i),n,N, j G No. 

Setting now r„+i = Ln+i,n,N and a^^N = Ln,n+d,N, the matrix element (n + d,n) of 
the equality fH3|) gives us relation fl3^ . 

The diagonal part of equation fHUj) now reads 



'r'm+i,N - ^li,N = 7^ + \^rn,N? " \am-d,N? + o(e ^^), Q < m < N - d. 
Adding up the first n + 1 of these equations gives us □ 
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In the special case ci = 1, we may set D = <C so that the error terms disappear. Then 
an^N = 2t2'"n+i,Ar and the recursion relation (JSlj) simplifies to 



which has the explicit solution r^^N = \J N(i-it^) ' ^ IN. 

Lemma 5.11. Let ((3'„)„6]No real polynomials with degg„ = n, go > and 
(55) zqn{z) = rn+iqn+i{z) + an-dqn-d{z), 

for some (i G N and some positive numbers and an, n & Nq. Then there uniquely exist 
real polynomials qj(^d+i)+r of degree j with 



qj{d+i)+r{z) = z''qj(^d+i)+r{z'^^^), <r <d, j G N, 



and these qj{^d+i)+r have exactly j different, real, positive zeros. 

Proof. Rewriting the recursion relation fl3^ in terms of the polynomials we get 

rn+iqn+i{z) = quiz) - an^dqn-d{z) for {d+ l)No + d and 

rn+iqn+i{z) = zqn{z) - an~dqn-d{z) for u E {d + l)No + d. 

Assume now by induction in j G Nq that all the zeros {^k,n)'lTi" of qn are positive numbers 
for j{d+l) <n< (j + l)(d + 1) and fulfil 

< 6,n+i < ^k+ijid+1) for j{d + 1) < n < j{d + 1) + d and I <k <j. 

(To simplify the notation, we set ^o,n = and ^degg„+i,n = oo for all n G Nq.) 
Then, for < k < j , 

sign(%-+i)(d+i)(^fc,j(d+i)+d)) = - sign(g(j+i)(d+i)(^fc+ij((i+i))), 

and so we have ^k,n < ^k+Uj+i){d+i) < ^k+i,n for all < A; < j and j{d + 1) < n < 
{j + l){d+l). 

Now we see recursively in 1 < r < that 

sign(%-+l)(d+l)+r(^A:,(j+l)(d+l)+r-l)) = - sign(g(j+i)(d+i)+^ (^fc ) ' < /c < j, 

verifying the induction hypothesis. □ 

Provided now that the error terms in recursion relation ( l52l) can be neglected (which is 
what we expect), we can apply Lemma [5 . 1 1 1 and see that in the limit N ^ oo the zeros of 
the orthonormal polynomials will lie on the tree-like graph {2; G C | z'^'^^ G [0,cx))}. As 
shown in Section 15.11 their density in the continuum limit can therefore be determined 
as the discontinuity of the Schwarz function of the polynomial curve parametrised by h. 
In the cases d = 1, where h parametrises an ellipse, and d = 2, where h parametrises a 
hypotrochoid, this can explicitly be done. 

2 2 

Proposition 5.12. The Schwarz function of the ellipse {{x,y) G | I3- + ^ = 1} with 



half-axis 61,62 > zs given by 
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Proof. Solving the quadratic equation 



4 



bj hi 

for 5(2;) and choosing the sign such that S{b\) = bi, leads us immediately to equation 

□ 

Proposition 5.13. Let C,^^ be the zeros of the orthonormal polynomial Qn^N for the 
potential 

V{z) = ^{\z\'-2t2iz^ + f)), 0<t2<l. 

to z 

Then, for t^ > small enough, each converging subsequence of the sequence (6c ) tends 

—n,N 

for N —>■ 00 and ^ ^ x to the measure a on H with da{y) = p{y) dy and 



2 



(57) piy) = j^^r^^bi -bl-y- ^[-xAF^.xAHli^^)- 

Here b\ and 62 denote the half-axis of the ellipse with vanishing first harmonic moment, 
second harmonic moment t2, and area irxtQ. 

Proof. Proposition 15.101 together with Lemma 15.111 tells us that the zeros all lie on the 
real axis. The density p is therefore given by 

_ 6S 

^ 27rixto ' 

where 6S is the discontinuity of the Schwarz function S of the ellipse, which was calculated 
in Proposition I5.12[ Here the branch cut of S has to be chosen on the part of the real 
axis where 6S is purely imaginary, that is on the interval [—^/bf — between 
the two foci of the ellipse. This directly brings us to equation (1571) . □ 

Proposition 5.14. 

(i) The function h{w) = rw + aw~'^, w G S'^, is for allr > and a G R with 2\a\ < r 
a parametrisation of a polynomial curve, the hypotrochoid. 

(ii) The Schwarz function S of this hypotrochoid is given by 

S{3rz) = (^{a^ - r^) 00 (^)' + {2a^ + r') 00 (^) + + 2r''^ ^, 
where 



2 ( , \ J 2 



uj{y) = ^i--[i + Vi-y) + ^i--[i-Vi-y 

and the branch cuts of the roots are chosen on the negative real axis. 

Proof. Because the derivative of h vanishes only for those w with rw^ = 2a, which for 
2\a\ < r all lie inside the unit disc, h is a biholomorphic map from the exterior of the unit 
disc and is therefore a parametrisation of a polynomial curve. 
Solving now h{w) = z ioi w leads us to the cubic equation 

~s ~ f ^ f z \^ a 
ir — ?ir)U! — 2(7 = D. -n = — I . n = \ — I . 



- 3pw - 2g = 0, p= j , q= 
for w = w — whose solution is of the form 



3r/ 2r 



w = SJ q + y'q'^ - p^ + y q - \/q'^ - p^. 
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Figure 1. Graph of ct(z) = p(f (i)^)^^)^ for z e [0,(f ar^)^] in the 



case to = 4, ^3 = 9 



Choosing the branch cuts of the cubic root such that the solution is continuous outside 
the unit disc, gives us the inverse function for h: 

So the Schwarz function S is given by 

2 2 

S{z) = + ah-\zf = -h~\zf + -h-\z) 



1 / ^ \2 



a \3r 



□ 



Proposition 5.15. Let be the zeros of the orthonormal polynomial qn^N for the 

potential 

v{z) = U\A'-^h{z' + z')), h>o. 

^0 

Then, provided the error terms in fl52l) are neglectable, for to > small enough, each 
converging subsequence of the sequence {6^ ) tends for N ^ oo and ^ —>■ x to the 
measure a. 



(58) 
where 

(59) 



1 



„ 1 - \ 1 - 72xtotl, a-- 
6t3 V V 3> ^2^3 



1 



1 - 72xtot^ 



Oh 5 



dz denotes the infinitesimal line element on the one- dimensional manifold {z E <C \ z E 

27 
4 



(0,f ar^)}, 



p{y) = {{u;4yy - u^{yf){T^ - a') - {u^{y) - u;4y)){r' + 2a')) 



27r(r2-2a2)' 



RANDOM NORMAL MATRICES AND POLYNOMIAL CURVES 



37 



and 

Proof. For small t^, the support of the equilibrium measure for the potential V is given 
by the polynomial curve 7 defined by the parametrisation flSOl) with d = 2. Solving now 
the biquadratic equation (ISTl) for r, we find, substituting xto for to, the relations ( !59l) . 

Proposition 15.141 provides us now with the Schwarz function S for the curve 7. The 
function 5* is discontinuous in those 2 G C, where 

l--l±yr^ <0, y = -(-) , 
y a V3r/ 

that is for y e (0, 1) resp. G (0, ^ar^). 

Because we know that all zeros of the orthonormal polynomials qn^N he on the set 
{z G C I z^ G [0,00)}, the limiting measure a is given through equation fHT]) leading us 
to formula ([58]). □ 
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